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OA , Abstract 

This paper considers the minimization of transmit power in Gaussian parallel interference channels, 
subject to a rate constraint for each user. To derive decentralized solutions that do not require 

o 

| equilibrium game. We obtain sufficient conditions that guarantee the existence and nonemptiness 

of the solution set to our problem. Then, to compute the solutions of the game, we propose two 
distributed algorithms based on the single user waterfilling solution: The sequential and the simul- 
taneous iterative waterfilling algorithms, wherein the users update their own strategies sequentially 
and simultaneously, respectively. We derive a unified set of sufficient conditions that guarantee the 
uniqueness of the solution and global convergence of both algorithms. Our results are applicable 
to all practical distributed multipoint-to-multipoint interference systems, either wired or wireless, 
where a quality of service in terms of information rate must be guaranteed for each link. 

Index Terms: Gaussian parallel interference channel, mutual information, game theory, general- 
ized Nash equilibrium, spectrum sharing, iterative waterfilling algorithm. 
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1 Introduction and Motivation 



The interference channel is a mathematical model relevant to many communication systems where 
multiple uncordinated links share a common communication medium, such as wireless ad-hoc networks 
or Digital Subscriber Lines (DSL). In this paper, we focus on the Gaussian parallel interference channel. 

A pragmatic approach that leads to an achievable region or inner bound of the capacity region 
is to restrict the system to operate as a set of independent units, i.e., not allowing multiuser encod- 
ing/decoding or the use of interference cancelation techniques. This achievable region is very relevant 
in practical systems with limitations on the decoder complexity and simplicity of the system. With 
this assumption, multiuser interference is treated as noise and the transmission strategy for each user is 
simply his power allocation. The system design reduces then to finding the optimum power distribution 
for each user over the parallel channels, according to a specified performance metric. 

Within this context, existing works [T]-|13j considered the maximization of the information rates 
of all the links, subject to transmit power and (possibly) mask constraints on each link. In [lj-|3j, 
a centralized approach based on duality theory [TU IE] was proposed to compute, under technical 
conditions, the largest achievable rate region of the system (i.e., the Pareto-optimal set of the achievable 
rates). In [4], sufficient conditions for the optimal spectrum sharing strategy maximizing the sum-rate 
to be frequency division multiple access (FDMA) were derived. However, the algorithms proposed 
in [I] -[4] are computationally expensive and cannot be implemented in a distributed way, require the 
full-knowledge of the system parameters, and are not guaranteed to converge to the global optimal 
solution. 

Therefore, in [5]-|13j. using a game-theory framework, the authors focused on distributed algorithms 
with no centralized control. In particular, the rate maximization problem was formulated as a strategic 
non-cooperative game, where every link is a player that competes against the others by choosing the 
transmission strategy that maximizes its own information rate [5j. Based on the celebrated notion 
of Nash Equilibrium (NE) in game theory |16j . an equilibrium for the whole system is reached when 
every player's reaction is "unilaterally optimal", i.e., when, given the rival players' current strategies, 
any change in a player's own strategy would result in a rate loss. In [6]-fT3], alternative sufficient 
conditions were derived that guarantee the uniqueness of the NE of the rate maximization game and 
the convergence of alternative distributed waterfilling based algorithms, either synchronous — sequential 
[6]-|ll| and simultaneous [T2] — or asynchronous [13J. 

The game theoretical formulation proposed in the cited papers, is a useful approach to devise totally 
distributed algorithms. However, due to possible asymmetries of the system and the inherent selfish 
nature of the optimization, the Nash equilibria of the rate maximization game in [S]-[T3] may lead 
to inefficient and unfair rate distributions among the links even when the game admits a unique NE. 
This unfairness is due to the fact that, without any additional constraint, the optimal power allocation 



2 



corresponding to a NE of the rate maximization game is often the one that assigns high rates to the 
users with the highest (equivalent) channels; which strongly penalizes all the other users. As many 
realistic communication systems require prescribed Quality of Service (QoS) guarantees in terms of 
achievable rate for each user, the system design based on the game theoretic formulation of the rate 
maximization might not be adequate. 

To overcome this problem, in this paper we introduce a new distributed system design, that takes ex- 
plicitly into account the rate constraints. More specifically, we propose a novel strategic non-cooperative 
game, where every link is a player that competes against the others by choosing the power allocation 
over the parallel channels that attains the desired information rate, with the minimum transmit power. 
We will refer to this new game as power minimization game. An equilibrium is achieved when every 
user realizes that, given the current power allocation of the others, any change in its own strategy 
would result in an increase in transmit power. This equilibrium is referred to as Generalized Nash 
Equilibrium (GNE) and the corresponding game is called Generalized Nash Equilibrium ProblemQ 

The game theoretical formulation proposed in this paper differs significantly from the rate maxi- 
mization games studied in [5]-[T3] . In fact, differently from these references, where the users are allowed 
to choose their own strategies independently from each other, in the power minimization game, the 
rate constraints induce a coupling among the players' admissible strategies, i.e., each player's strategy 
set depends on the current strategies of all the other players. This coupling makes the study of the 
proposed game much harder than that of the rate maximization game and no previous result in [5]-|13j 
can be used. Recently, the calculation of generalized Nash equilibria has been the subject of a renewed 
attention also in the mathematical programming community, see for example [17j-|21j. Nevertheless, 
in spite of several interesting advances [21], none of the game results in the literature are applicable to 
the power minimization game. 

The main contributions of the paper are the following. We provide sufficient conditions for the 

nonemptiness and boundedness of the solution set of the generalized Nash problem. Interestingly, 

these sufficient conditions suggest a simple admission control procedure to guarantee the feasibility 

of a given rate profile of the users. Indeed, our existence proof uses an advanced degree-theoretic 

result for a nonlinear complementarity problem in order to handle the unboundedness of the users' 

rate constraints. We also derive conditions for the uniqueness of the GNE. Interestingly, our sufficient 

conditions become also necessary in the case of one subchannel. To compute the generalized Nash 

solutions, we propose two alternative totally distributed algorithms based on the single user waterfilling 

solution: The sequential IWFA and the simultaneous IWFA. The sequential IWFA is an instance of the 

Gauss-Seidel scheme: The users update their own strategy sequentially, one after the other, according 

1 According to recent use, we term generalized Nash equilibrium problem a Nash game where the feasible sets of the 
players depend on the other players' strategy. Such kind of games have been called in various different ways in the 
literature, for example social equilibrium problems or just Nash equilibrium problems. 
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to the single user waterfilling solution and treating the interference generated by the others as additive 
noise. The simultaneous IWFA is based on the Jacobi scheme: The users choose their own power 
allocation simultaneously, still using the single user waterfilling solution. Interestingly, even though the 
rate constraints induce a coupling among the feasible strategies of all the users, both algorithms are still 
totally distributed. In fact, each user, to compute the waterfilling solution, only needs to measure the 
power of the noise plus the interference generated by the other users over each subchannel. It turns out 
that the conditions for the uniqueness of the GNE are sufficient for the convergence of both algorithms. 
Our convergence analysis is based on a nonlinear transformation that turns the generalized game in 
the power variables into a standard game in the rate variables. Overall, this paper offers two major 
contributions to the literature of game-theoretic approaches to multiuser communication systems: (i) 
a new noncooperative game model is introduced for the first time that directly addresses the issue of 
QoS in such systems, and (ii) a new line of analysis is introduced in the literature of distributed power 
allocation that is expected to be broadly applicable for other game models. 

The paper is organized as follows. Section [2] gives the system model and formulates the power 
minimization problem as a strategic non-cooperative game. Section [3] provides the sufficient conditions 
for the existence and uniqueness of a GNE of the power minimization game. Section [4] contains the 
description of the distributed algorithms along with their convergence conditions. Finally, Section [5] 
draws the conclusions. Proofs of the results are given in the Appendices IAHfI 

2 System Model and Problem Formulation 

In this section we clarify the assumptions and the constraints underlying the system model and we 
formulate the optimization problem explicitly. 

2.1 System model 

We consider a Q-user Gaussian iV-parallel interference channel. In this model, there are Q transmitter- 
receiver pairs, where each transmitter wants to communicate with its corresponding receiver over a 
set of N parallel subchannels. These subchannels can model either frequency-selective or flat-fading 
time-selective channels [9]. Since our goal is to find distributed algorithms that do not require neither 
a centralized control nor a coordination among the links, we focus on transmission techniques where no 
interference cancelation is performed and multiuser interference is treated as additive colored noise from 
each receiver. Moreover, we assume perfect channel state information at both transmitter and receiver 
side of each link(^| each receiver is also assumed to measure with no errors the power of the noise plus 
the overall interference generated by the other users over the N subchannels. For each transmitter q, 
2 Note that each user is only required to known its own channel, but not the channels of the other users. 
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the total average transmit power over the N subchannels is (in units of energy per transmitted symbol) 

1 N 

k=l 

where p q (k) denotes the power allocated by user q over the subchannel k. 

Under these assumptions, invoking the capacity expression for the single user Gaussian channel — 
achievable using random Gaussian codes from all the users — the maximum information rate on link q 
for a specific power allocation is |23 J 



JV 

^(Pg,P-«) = X)log(l + sinr ff (fc)), (2) 

k=l 

with smr q (k) denoting the Signal-to-Interference plus Noise Ratio (SINR) of link q on the fc-th sub- 
channel: 

sinr (k) 4 \H qq (k)\ 2 p q (k) 

11 V 2 q (k) + Zr^ q \H q r(k)\ 2 Pr(ky 

where \H qr (k)\ 2 is the power gain of the channel between destination q and source r; cr 2 (k) is the 
variance of Gaussian zero mean noise on subchannel k of receiver q; and p q = {p q (k))^ =1 is the power 
allocation strategy of user q across the N subchannel, whereas p_ g = (p T ) r ^ q contains the strategies 
of all the other users. 

2.2 Game theoretic formulation 

We formulate the system design within the framework of game theory [25j [26], using as desirability 
criterion the concept of GNE, see for example [16j [27j. Specifically, we consider a strategic non- 
cooperative game, in which the players are the links and the payoff functions are the transmit powers 
of the users: Each player competes against the others by choosing the power allocation (i.e., its strategy) 
that minimizes its own transmit power, given a constraint on the minimum achievable information rate 
on the link. A GNE of the game is reached when each user, given the strategy profile of the others, does 
not get any power decrease by unilaterally changing its own strategy, still keeping the rate constraint 
satisfied. Stated in mathematical terms, the game has the following structure: 

sr = {n,{& q ( P „ q )} qen ,{p q } qeQ } , (4) 



where $7 = {1,2, . . . , Q} denotes the set of the active links, & q (p- q ) C is the set of admissible 
power allocation strategies p q G HP q (p- q ) of user q over the subchannels M = {1, . . . , N}, defined as 

&» q (p- q ) = {x 9 G : R q (x q ,p- q ) > R*} . (5) 



3 Observe that a GNE is obtained if each user transmits using Gaussian signaling, with a proper power allocation. 
However, generalized Nash equilibria achievable using non-Gaussian codes may exist. In this paper, we focus only on 
transmissions using Gaussian codebooks. 
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with Rq(p q ,p_ q ) given in j2]), and R q denotes the minimum transmission rate required by each user, 
which we assume positive without loss of generality. In the sequel we will make reference to the vector 
R* = (R q ) q=1 as to the rate profile. The payoff function of the q-th player is its own transmit power 
P q , given in ([TJ). Observe that, because of the rate constraints, the set of feasible strategies ^ g q (p- q ) 
of each player q depends on the power allocations p- q of all the other users. 

The optimal strategy for the q-th player, given the power allocation of the others, is then the 
solution to the following minimization problem 

N 

minimize > p a (k) 

ti . (6) 

subject to p q G &> q (p- q ) 

where 3? q {p- q ) is given in ((5|). Note that, for each q, the minimum in j6]) is taken over p g , for a fixed 
but arbitrary P- q . Interestingly, given P- q , the solution of (JSJ) can be obtained in "closed" form via 
the solution of a singly-constrained optimization problem; see [28] for an algorithm to implement this 
solution in practice. 

Lemma 1 For any fixed and nonnegative p~ q , the optimal solution p q = {p q (k)}^ =1 of the optimization 
problem fSjj exists and is unique. Furthermore, 

Pq = WF 9 (Pl» • • • I P«-l» Pff+li • • • , Pq) = WF,(p_ 9 ) , (7) 
where the waterfilling operator WF g (•) is defined as 

( - \H qr (k)\ 2 p r (k)\ + 

[WF 9 (p- q )] k = J X q - 2 , k € Af, (8) 

y \ H qq{k)\ J 

with (x) + = max(0,x) and the water-level X q chosen to satisfy the rate constraint R q (p*,p- q ) = R q , 
with Rq(Pq,P- q ) given in |^). 

The solutions of the game Sf in if they exist, are the Generalized Nash Equilibria, formally 
defined as follows. 



Definition 2 A feasible strategy profile p* = (p*)® =1 is a GNE of the game if 

N N 

< 5>,(fc), Vp, G & q (p%), Vq G 0. (9) 

k=l k=l 

According to Lemma [U all the Generalized Nash Equilibria of the game must satisfy the condition 
expressed by the following Corollary. 
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Corollary 3 A feasible strategy profile p* = (Pg)^ =1 is a GNE of the game 5f if and only if it satisfies 
the following system of nonlinear equations 



Pg = WF 9(Pl'---' P 9-l' P 3+l' 



Vg G n, 



(10) 



with WF g (•) defined in 



Given the nonlinear system of equations (|10p . the fundamental questions we want an answer to are: 
i) Does a solution exist, for any given users' rate profile? ii) // a solution exists, is it unique? iii) How 
can such a solution be reached in a totally distributed way? 

An answer to the above questions is given in the forthcoming sections. 

3 Existence and Uniqueness of a Generalized Nash Equilibrium 

In this section we first provide sufficient conditions for the existence of a nonempty and bounded 
solution set of the Nash equilibrium problem (Jl|). Then, we focus on the uniqueness of the equilibrium. 



3.1 Existence of a generalized Nash equilibrium 

Given the rate profile R* = (Rq)? =1 , define, for each k G J\f, the matrix Z/-(R*) G M. Q 

\H n (k)\ 2 -(e R * - 1) \H 12 (k)\ 2 ■ ■ ■ -(e R * - 1) \H 1Q (k)\ 2 



as 



Z fc (R*) 



(e*S -1) \H 21 (k)f 



I #22 (k) | 



-(e R Q - 1) \H Q1 (k)\ z -{e K Q - 1) |F Q2 (A;)| 
We also need the definition of P-matrix, as given next. 



(e*2_l) \H 2Q {k)Y 



\H QQ {k)Y 



fill 



Definition 4 A matrix A G M. NxN is called Z-matrix if its off-diagonal entries are all non- positive. 
A matrix A G M. NxN is called P-matrix if every principal minor of A is positive. 

Many equivalent characterizations for a P-matrix can be given. The interested reader is referred to 
pUJ E2] for more details. Here we note only that any positive definite matrix is a P-matrix, but the 
reverse does not hold. 

Sufficient conditions for the nonemptiness of a bounded solution set for the game 5f are given in 
the following theorem. 

Theorem 5 The game with rate profile R* = (R*)® =1 > admits a nonempty and bounded solution 
set i/Zfc(R*) is a P-matrix, for all k G N, with Zfc(R*) defined in / TO]) . Moreover, any GNE p* = 



7 



(jPq)9 = i is such that 



( Plik) \ 



I Pi(k) \ 



< 



( a 2 {k)(e R i - 1) \ 



, k G A/\ 



:i2) 



\ a 2 Q {k) (e R Q - 1) J 



(Z^R^))- 1 

Proof. See Appendix lAl ■ 
A more general (but less easy to check) result on the existence of a bounded solution set for the 

game is given by Theorem [23] in Appendix [Al 

We now provide alternative sufficient conditions for Theorem [5] in terms of a single matrix. To this 

end, we first introduce the following matrix 

-(e R * - l)/?^ ax • • • -(e R * - 1) /3^ ax 



Z max (R*) — 



-(e R z - l)f3. 



max 
'21 



(e R 2 - l)/3^ ax 



where 



max 
Q2 



(13) 



qr 



max 



keN \H„{k)Y 



Vr q , q £ Q. 



(14) 



We also denote by e* 1 — 1 the Q-vector with g-th component e R i — 1, for q = 1, . . . , Q. Then, we have 
the following corollary 

Corollary 6 // Z max (R*) in [Hfy is a P-matrix, then all the matrices {Zfc(R*)} defined in {77p are 
P-matrices. Moreover, any GNE p* = (p*)^_ 1 of the game 5f satisfies 



P*q(k) <pJk) 



max<7^(/c) 



H m {k)Y 



(Z max (R^))~ 1 ( e R ~ - 1 ) , \/q G f2, Vfc G M. (15) 



Proof. See Appendix iBl ■ 
To give additional insight into the physical interpretation of the existence conditions of a GNE, 
we make explicit the dependence of each channel (power) gain \H qr (k)\ 2 on its own source-destination 
distance d qr by introducing the normalized channel gain \H qr (k)\ 2 = \H qr (k)\ 2 dq r , where 7 is the path 
loss exponent. We have the following corollary. 



Corollary 7 Sufficient conditions for the matrices {Z/%(R*)} defined in ^11}) to be P-matrices are: 



1 7? 



H qr (k)\ 2 d, 



2 Al 
qq 



< 



Vr £Q, \fk G M. 



(16) 
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Proof. The proof comes directly from the sufficiency of the diagonally dominance property [321 Defi- 
nition 2.2.19] for the matrices Z fc (R*) in ((TTJl to be P-matrices [3H Theorem 6.2.3] ■ 



Remark 8 A physical interpretation of the conditions in Theorem [5] (or Corollary [7J is the following. 
Given the set of channels and the rate constraints, a GNE of is guaranteed to exist if multiuser 
interference is "sufficiently small" (e.g., the links are sufficiently far apart). In fact, from (fT6j) . which 
quantifies the concept of small interference, one infers that, for any fixed set of (normalized) channels 
and rate constraints, there exists a minimum distance beyond which an equilibrium exists, correspond- 
ing to the maximum level of interference that may be tolerated from each user. The amount of such 
a tolerable multiuser interference depends on the rate constraints: the larger the required rate from 
each user, the lower the level of interference guaranteeing the existence of a solution. The reason why 
an equilibrium of the game £f might not exist for any given set of channels and rate constraints, is 
that the multiuser system we consider is interference limited, and thus not every QoS requirement is 
guaranteed to be feasible. In fact, in the game , each user acts to increase the transmit power to 
satisfy its own rate constraint; which leads to an increase of the interference against the users. It turns 
out that, increasing the transmit power of all the users does not guarantee that an equilibrium could 
exist for any given rate profile. 

Observe that conditions in Theorem [5] also provide a simple admission control procedure to check 
if a set of rate constraints is feasible: under these conditions indeed, one can always find a finite power 
budget for all the users such that there exists a GNE where all the rate constraints are satisfied. 

3.2 Uniqueness of the Generalized Nash Equilibrium 

Before providing conditions for the uniqueness of the GNE of the game , we introduce the following 
intermediate definitions. For any given rate profile R* = (R*)® =1 > 0, let B(R*) G m9 x ® be defined 



as 




if q = r, 



otherwise 



(17) 



where 



qr 



max 



/ \H qr (k)\ 2 a 2 Ak) + Er Vr \ H rr ' I ' Pr > (g 

keM \ \H rr (k)\ 2 (7|(fe) 



) 



(18) 



with p r /(k) defined in (fT2l) 



We also introduce x an d p, defined respectively as 




(19) 



with /3™ ax given in (fHj) . and 



A 

P = 




- 1 



- 1 



(20) 
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with 

^max-max-Rq, and R^ in = mmR* (21) 
Sufficient conditions for the uniqueness of the GNE of the game are given in the following theorem. 

Theorem 9 Given the game <S with a rate profile R* = (R*)® =1 > 0, assume that the conditions of 
Theorem^ are satisfied. If, in addition, B(R*) in SI 7\) is a P-matrix, then the game admits a unique 
GNE. 

Proof. See Appendix iDl ■ 
More stringent but more intuitive conditions for the uniqueness of the GNE are given in the following 
corollary. 

Corollary 10 Given the game W with rate profile R* = (R*)® =1 > 0, assume that 

< X < 1, (22) 

so that a GNE for the game ^ is guaranteed to exist, with x defined in / flP)) . Then, the GNE is unique 
if the following conditions hold true 

V/3- ax { ( max + [max ( max 

with p defined in IW^ . 

In particular, when cr^(k) = <r?(ra), Vr, q G and V/c, n E jV, conditions i23\) become 

(\H qr (k)\ 2 ) dr? r 7 w „ 

gs{fewp}< < ?fr? v ' £ "' (24) 



^- 1 )} < ^' V96fi ' (23) 



with 



7 = -r* : <!• (25 

e max — 1 _ _9R*i 

+ max { e s — e « j 



max {e R *i — e 2i? «} 



e R min — 1 g eO 

Proof. See Appendix [El ■ 
3.3 On the conditions for existence and uniqueness of the GNE 

It is natural to ask whether the sufficient conditions as given by Theorem [5] (or the more general ones 
given by Theorem [23] in Appendix A) are tight. In the next proposition, we show that these conditions 
become indeed necessary in the special case of N = 1 subchannel. 

Proposition 11 Given the rate profile R* = (R*) q= ^, the following statements are equivalent for the 
game *S when N = lQ 



4 In the case of N = 1, the power allocation p q (k) — p q of each user, the channel gains \H rq {k)\ = \H rq \ and the noise 
variances <J q (k) = a 2 are independent on index k. Matrix Zfc(R*) = Z(R*) is defined as in l|lip . where each \H rq (k)\ 2 is 
replaced by |-HV 9 | 2 . 
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(a) The problem has a solution for some (all) (o"g)^ =1 > 0. 

(b) The matrix Z(R*) is a P-matrix. 

If any one of the above two statements holds, then the game 'S has a unique solution that is the unique 
solution to the system of linear equations: 

\H qq f P q -(e R « \H qr \ 2 Pr = a 2 g (e R «-l) \/q G Q. (26) 

r+q 

Proof. See Appendix ■ 

Remark 12 Proposition [TT] also shows that it is, in general, very hard to obtain improved sufficient 
condition for the existence and boundedness of solutions to the problem with N > 1, as any such 
condition must be implied by condition (b) above for the the 1-subchannel case, which, as shown by 
the proposition, is necessary for the said existence, and also for the uniqueness as it turns out. 

Remark 13 Observe that, when N = 1, the game <S leads to classical SINR based scalar power control 
problems in flat-fading CDMA (or TDMA/FDMA) systems, where the goal of each user is to reach a 
prescribed SINR (see with the minimum transmit power P q [22] • In this case, given the rate profile 
R* = (R*)y =1 and N = 1, the SINR target profile sinr* = (sinr*)^ =1 , as required in classical power 
control problems [29J, can be equivalently written in terms of R*as 

sinr* = e~ R « - 1, q G U, (27) 

and the Nash equilibria p* = (pg)® = i of the game §f become the solutions of the following system of 
linear equations 

' ufsinr* ^ 



Z(R*)p* 



(2E 



V °Q s ' nr Q / 



Interestingly, the necessary and sufficient condition (b) given in Proposition [TT] is equivalent to that 
known in the literature for the existence and uniqueness of the solution of the classical SINR based 
power control problem (see, e.g., [29J). Moreover, observe that, in the case of N = 1, the solution of 
the game £f, coincides with the upper bound in (fT2l) . 

Numerical example. Since the existence and uniqueness conditions of the GNE given so far depend 
on the channel power gains {\H qr (k)\ 2 }, there is a nonzero probability that they are not satisfied 
for a given channel realization drawn from a given probability space and rate profile. To quantify the 
adequacy of our conditions, we tested them over a set of channel impulse responses generated as vectors 
composed of L = 6 i.i.d. complex Gaussian random variables with zero mean and variance equal to the 
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square distance between the associated transmitter-receiver links (multipath Rayleigh fading model). 
Each user transmits over a set of N = 32 subcarriers. We consider a multicell cellular network as 
depicted in Figure [TJi) , composed of 7 (regular) hexagonal cells, sharing the same band. Hence, the 
transmissions from different cells typically interfere with each other. For the simplicity of representation, 
we assume that in each cell there is only one active link, corresponding to the transmission from the 
base station (BS) to a mobile terminal (MT). According to this geometry, each MT receives a useful 
signal that is comparable, in average sense, with the interference signal transmitted by the BSs of 
two adjacent cells. The overall network can be modeled as a 7-users interference channel, composed 
of 32 parallel subchannels. In Figure [lb), we plot the probability that existence (red line curves) and 
uniqueness (blue line curves) conditions as given in Theorem [5] and Theorem [U respectively, are satisfied 
versus the (normalized) distance d € [0,1) [see Figure Hk)], between each MT and his BS (assumed 
to be equal for all the MT/BS pairs). We considered two different rate profiles, namely R* = 1 
bit /symb/ sub channel (square markers) and R* = 2 bit/symb/subchannel (cross markers), Vg G fl. As 
expected, the probability of existence and uniqueness of the GNE increases as each MT approaches his 
BS (i.e., d — ► 1), corresponding to a decrease of the intercell interference. 



4 Distributed Algorithms 

The game was shown to admit a GNE, under some technical conditions, where each user attains 
the desired information rate with the minimum transmit power, given the power allocations at the 
equilibrium of the others. In this section, we focus on algorithms to compute these solutions. Since we 
are interested in a decentralized implementation, where no signaling among different users is allowed, 
we consider totally distributed algorithms, where each user acts independently to optimize its own 
power allocation while perceiving the other users as interference. More specifically, we propose two 
alternative totally distributed algorithms based on the waterfilling solution in J7|), and provide a unified 
set of convergence conditions for both algorithms. 

4.1 Sequential iterative waterfilling algorithm 

The sequential Iterative Waterfilling Algorithm (IWFA) we propose is an instance of the Gauss-Seidel 
scheme (by which, each user's power is sequentially updated [22]) based on the mapping (J7]): Each 
player, sequentially and according to a fixed updating order, solves problem j6|), performing the single- 
user waterfilling solution in J7]). The sequential IWFA is described in Algorithm 1. 
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Algorithm 1: Sequential Iterative Waterfilling Algorithm 



Set p q 



(0) 



any nonnegative vector; 



f or n = : Number of iterations 



(n+l) 



= < 



WF, 



q (p-g) , if (» + 1) mod Q = q 



Vg G ft; 



(29) 




otherwise 



end 



The convergence of the algorithm is guaranteed under the following sufficient conditions. 

Theorem 14 Assuming Number_ of_ iterations = oo, the sequential IWFA, described in Algorithm 
1, converges linearly^ to the unique GNE of the game & , if the conditions of Theorem^ are satisfied. 



Remark 15 Observe that the convergence of the algorithm is guaranteed under the same conditions 
obtained for the uniqueness of the solution of the game. As expected, the convergence is ensured if the 
level of interference in the network is not too high. 

Remark 16 The main features of the proposed algorithm are its low- complexity and distributed na- 
ture. In fact, despite the coupling among the users' admissible strategies due to the rate constraints, the 
algorithm can be implemented in a totally distributed way, since each user, to compute the waterfilling 
solution j7]), only needs to locally measure the interference-plus-noise power over the N subchannels 
[see j3])] and waterfill over this level. 

Remark 17 Despite its appealing properties, the sequential IWFA described in Algorithm 1 may suffer 

from slow convergence if the number of users in the network is large, as we will also show numerically 

in Section 14.21 This drawback is due to the sequential schedule in the users' updates, wherein each 

user, to choose its own strategy, is forced to wait for all the other users scheduled before it. It turns 

out that the sequential schedule, as in Algorithm 1, does not really gain from the distributed nature 

of the multiuser system, where each user, in principle, is able to change its own strategy, irrespective 

of the update times of the other users. Moreover, to be performed, the sequential update requires a 

centralized synchronization mechanism that determines the order and the update times of the users. 

We address more precisely this issue in the next section. 

6 A sequence {x n } is said to converge linearly to x* if there is a constant < c < 1 such that ||2;n+i — x*\\ < c\\x n — x*\\ 
for all n > n and some n £ N. 



Proof. See Appendix [H 
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4.2 Simultaneous iterative waterfilling algorithm 

To overcome the drawback of the possible slow speed of convergence, we consider in this section the 
simultaneous version of the IWFA, called simultaneous Iterative-waternlling Algorithm. The algorithm 
is an instance of the Jacobi scheme [22j: At each iteration, the users update their own PSD simulta- 
neously, performing the waterfilling solution ([3), given the interference generated by the other users in 
the previous iteration. The simultaneous IWFA is described in Algorithm 2. 



Algorithm 2: Simultaneous Iterative Waterfilling Algorithm 

Set Pq = any nonnegative vector, Vg € fi; 
for n = : Number_of_ iterations 

p(™ +1) = WF 9 (pl n) , . . .,P q %P q % • • - , P W) , \fq G Q, (30) 

end 



Interestingly, (sufficient) conditions for the convergence of the simultaneous IWFA are the same as 
those required by the sequential IWFA, as given in the following. 

Theorem 18 Assuming Number_ of_ iterations = oo, the simultaneous IWFA, described in Algo- 
rithm 2, converges linearly to the unique GNE of the game ^ , if the conditions of Theorem are 
satisfied. 

Proof. See Appendix [Fl ■ 

Remark 19 Since the simultaneous IWFA is still based on the waterfilling solution j7]), it keeps the 
most appealing features of the sequential IWFA, namely its low- complexity and distributed nature. In 
addition, thanks to the Jacobi-based update, all the users are allowed to choose their optimal power 
allocation simultaneously. Hence, the simultaneous IWFA is expected to be faster than the sequential 
IWFA, especially if the number of active users in the network is large. 

Numerical Example. As an example, in FigureEl we compare the performance of the sequential and 
simultaneous IWFA, in terms of convergence speed. We consider a network composed of 10 links and 
we show the rate evolution of three of the links corresponding to the sequential IWFA and simultaneous 
IWFA as a function of the iteration index n as defined in Algorithms 1 and 2. In Figured) we consider 
a rate profile for the users with two different classes of service; whereas in Figure [2b) the same target 
rate for all the users is required. As expected, the sequential IWFA is slower than the simultaneous 
IWFA, especially if the number of active links Q is large, since each user is forced to wait for all the 
other users scheduled before updating its power allocation. 
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5 Conclusions 



In this paper we have considered the distributed power allocation in Gaussian parallel interference 
channels, subject to QoS constraints. More specifically, we have proposed a new game theoretic for- 
mulation of the power control problem, where each user aims at minimizing the transmit power while 
guaranteeing a prescribed information rate. We have provided sufficient conditions for the nonempti- 
ness and the boundedness of the solution set of the Nash problem. These conditions suggest a simple 
admission control procedure to check the feasibility of any given users' rate profile. As expected, there 
exists a trade-off between the performance achievable from each user (i.e., the achievable information 
rate) and the maximum level of interference that may be tolerated in the network. Under some ad- 
ditional conditions we have shown that the solution of the generalized Nash problem is unique and 
we have proved the convergence of two distributed algorithms: The sequential and the simultaneous 
I WFAs. Interestingly, although the rate constraints induce a coupling among the feasible strategies of 
the users, both algorithms are totally distributed, since each user, to compute the waterfilling solution, 
only needs to measure the noise-plus-interference power across the subchannels. Our results are thus 
appealing in all the practical distributed multipoint-to-multipoint systems, either wired or wireless, 
where centralized power control techniques are not allowed and QoS in terms of information rate must 
be guaranteed for each link. 

One interesting direction that is worth of further investigations is the generalization of the proposed 
algorithms to the case of asynchronous transmission and totally asynchronous updates among the users, 
as did in [13] for the rate maximization game. 

6 Appendices 

A Proof of Theorem [5] 

We derive Theorem [5] as a corollary to the more general Theorem [23] below. In order to prove this 
theorem we need several preliminary concepts and results though, as given next. 

A.l Noiseless game 

We rewrite first the KKT optimality conditions of the Nash problem ((6)) as a Mixed nonlinear Comple- 
mentarity Problem (MNCP) [EJUS]- Denoting by [i q ^ the multipliers of the nonnegativity constraints 
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and by X q the multipliers of the rate constraint, the KKT conditions of problems (JSJ) can be written as: 



\H (k)\ 2 

1 - tie k ~ \ n o = 0> VA; G M, Vq G ft, 

' ^k) + E Q r=1 \H qr (k)\ 2 Mk) 



< p q (k) _L ng }k > 0, V/c G AT, G O, 

o<a 9 i Viog|i+ \H qq (k)\ 2 Pq (k) )-m >0 , v q en, 

(31) 

where a _L b means the two scalars (or vectors) a and b are orthogonal. Observe that each X q > 0; 
otherwise complementarity yields p q (k) = for all k G M, which contradicts the rate constraints. 

Eliminating the multipliers {/J. q ^} corresponding to the nonnegativity constraints and making some 
obvious scaling, the KKT conditions in (|3T1) are equivalent to the following MNCP: 

Q 

< p q {k) ± a 2 q (k) + \H qr {k)\ 2 p r (k) - \H qq (k)\ 2 \ q > 0, Vk G AT, Vg G fi. 

(32) 




To proceed further we introduce an additional game, which has the same structure of the game , 
except for the players' payoff functions, defined as in PJ, but with cr 2 (k) = for all k G M and g 6 !1. 
We will refer to this game as the noiseless game. Although the noiseless game does not correspond to 
any realistic communication system, it will be shown to be instrumental in understanding the behavior 
of the original game Sf when all cr 2 (k) > 0. 

Note that the conditions cr 2 (k) > ensure that all the users' rates R q (p q , P- 9 ) in J2)) of the <S in j6]) 
are well-defined for all nonnegative p = (p q )9 = i, with p q = {p q (k))^ =1 . Nevertheless, when o~ 2 {k) = 0, 



the players' payoff functions j R q (p q , p~ q ) of the noiseless game still remain well-defined as long as 
Yl?=iPr(k) > 0, provided that we allow for a rate equal to oo. Most importantly, the MNCP (|32j) is 
well defined for all nonnegative cr 2 (k), including the case when a 2 (k) = for all k G N and q G The 
latter observation motivates the following definition. 

Definition 20 A set of user powers p = (p q ) q= i, with p q = (p q (k))^ =1 , is said to be an almost GNE 
of the noiseless game if there exists a set of nonnegative scalars {v q }® =1 such that 

Q 

< Pq (k) ± \H qr (k)\ 2 p r (k) - \H qq (k)\ 2 v q > 0, \/k G M, Vg G Q. (33) 

r=l 

We call these solutions noiseless almost equilibria, and denote the set of noiseless equilibria by N£q. □ 

The set N£$ of users' noiseless almost equilibria constitutes a closed, albeit not necessarily convex, 
cone in the space of all users' powers. A noteworthy point about such an almost equilibrium is the 



6 With a slight abuse of notation, we use the same symbol to denote the payoff functions of the players in the game C S 
in J6J and in the noiseless game. 
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following simple property, which asserts that in the noiseless game, every subchannel will be used by 
at least one user. 

Proposition 21 If p G M£q, and p / 0, then Yl?=iPr(k) > for all keJ\f. 

Proof. Let p G M£o be such that p qo (ko) > for some pair (qo,ko). By complementarity (see j33j)), 
we have 

Q 

\H qoqo (k )\ 2 v qo = ^2\H qor (k )\ 2 Pr (k ) > 0, (34) 
r=l 

which implies 

Q 

J2\H qor (k)\ 2 Pr (k) > \H qoqo (k)\ 2 v qo > 0, Mk G J\f, (35) 

r=l 

since the \H qr (k)\ 2 are all positive. Equation ([35]) clearly implies: 

Q 

Y,Pr(k) > 0, \/k G Af, 
r=l 

as claimed by the proposition. ■ 
A. 2 The noiseless asymptotic cone 

Another mathematical concept we need is that of asymptotic direction of a (nonconvex) set that we 
borrow from recession analysis |30j . 

Given the game with rate profile R* = (R q )® =1 , it is possible that the sets of powers p that allow 
the user to achieve this rate be unbounded. In essence, the asymptotic consideration below aims at 
identifying such unbounded user powers. Specifically, we consider the following noisy nonconvex level 
set of users' powers corresponding to R*: 

= { P > o = £ >°s (i + ^fCTM^) ) =W»\ <*> 

where a = (<r 9 )^_ 1 , with a q = {o~ 2 (k))^ =1 . The asymptotic cone of 7 7<T (R*), denoted by V^(R*), is 
the cone (not necessarily convex) of directions d = ( d q )g =1 , with d q = ( d q (k))^ =1 , such that 

d = lim — (37) 

for some sequence of scalars { t u }^ =1 tending to oo and some sequence of powers { p v }^ = i such that 
p u G P a {W) for all v. 

It is known [30l Proposition 2.1.2] that P cr (R*) is bounded if and only if P£(R*) = {0}. Whereas 
the individual sets 'P <T (R*) are dependent on o- 2 {k), it turns out that the asymptotic cones of all these 
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sets are the same and equal to the following noiseless level set of users' powers: 



P U (R*) = { 



fe:Er=l^( & )>° 

where by convention the vacuous summation is defined to be zero (i.e., by definition, 'P°(R*) contains 
the origin). The claim about the equality of V^(R ) for a fixed R is formally stated and proved in 
the result below. 

Proposition 22 For any a > 0, P£(R*) = V°(R*). 

Proof. Let er > be arbitrary. Let d G V£>(R*) and {(p u ,T u )} be a sequence satisfying the definition 
of d in ([371) . For each q G £1 and all u, we have 

/ \ 



N 



R* > ^log 



fc=l 



1 + 



\H q Jk)\ 2 p»Jk) 



a 2 (k) +J2\H qr m 2 P?(k) 



( 



fc:£?=idr(*)>0 



\ 



1 + 



|ff w (*0W*0 



(38) 



fc:E? =1 d r (fc)>0 



H qq (k)\ 



2 



1 + 



a 2 q (k) 



V 



J 



Taking the limit v — > oo establishes the inclusion P£(R*) C P°(R*). 

Conversely, it is clear that G 7 3<T (R*). Let d be a nonzero vector in 7 ?0 (R*). For any scalar 
> 0, we have, for all q G fi, 



5>g 1 + 



0|F™(AO| 2 d,(fc) 



fe=l 



V 2 q (k) + 6Er H \Hqr(kWd r (k) J 

fc:E?=irfr(fc)>0 ^ 



e\H qq (k)\ 2 d q (k) \ 

^ 2 q (k) + 9E^ q \H q r(k)\ 2 d r (k) ) (39) 
\H qq (k)\ 2 d q (k) \ 



' £^j 0g V + ^q\Hqr(k)\ 2 d r (k) 



Hence #d G P°"(R*) for all 9 > 0. It follows that d G P£(R*) 
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We are now ready to introduce the key object in our proof of Theorem [23], the cone MSq fl 7- ,0 (R*), 
which by Proposition l22l is equal to {0} U AA£o(R*), where 



A/* (R*) 4 J d G A^o \ {0} : f log f i + = ^#^#77, ) < ^ Vg € n 



fc=l 



£ r ^|iM*OI 2 <*r(fc) 



(40) 



= A/"£ n P°(R*) \ {0}. 
Notice that d G AA£ (R*) implies £ r -^ l#gr(k)| 2 d r (£;) > for all q £ O, and all k G jV. 

A. 3 Existence results 

With the above preparation, we are ready to present our main existence theorem. The emptiness of 
the set MS o(R*) defined in ([4*01) turns out to provide a sufficient condition for the MNCP ([32]) to have 
a nonempty bounded solution set. 

Theorem 23 Given the game <S with rate profile R* = (R*)^ =1 > 0, if M£o(R*) = %, then the game 
has a nonempty and bounded solution set, for all a > 0. 

Proof. We first note that the KKT conditions of the Nash problem defined in ([6|) are equivalent to 
the following nonlinear complementarity problem (NCP) (see (|31|) and comments thereafter) 

Q 

0< Pq (k) J_ a 2 q (k) + Y,\H qr (k)\ 2 p r (k)-\H qq (k)\ 2 X q > 0, Vk G M, G O, 

r=l 

(41) 

In turn, to show that (I4*T1) has a solution, it suffices to prove that the solutions of the augmented NCP 

Q 

0< Pq (k) ± a 2 q {k) + Y,\H qr {k)\ 2 p r {k)- \H gq (k)\ 2 X q + r Pq (k) > 0, Vk G AA, G fi, 

r=l 



< A 9 1 ^bf 



/ \ 

1+ \H qq (k)\ 2 Pq (k) 



Rq + t \ q > o, v<7 g n, 



a 2 (k) + Y,\H qr (k)\ 2 p r (k) 

\ r^q J 

for all t > are bounded [17j Theorem 2.6.1]. 

We show the latter boundedness property by contradiction. Assume that for some sequence of 
positive scalars {t u }, a sequence of solutions {(p u ,r u )} exists such that each pair (p u ,r u ) satisfies: 

Q 

Q<P v q {k) ± a 2 q (k) + Y, \H qr (k)\ 2 p»{k) - \H qq (k)\ 2 y q + r u p-{k) > \/keM,Vqen, 

r=l 

(42) 

N ( \H (k)\ 2 v u (k) \ 

0SA « x ^ log ( 1+ WCT^ra l- fl ' +T - A '- ' v ^"- 



fc=i 
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and that 



lim 

v— >oo 



IP 



q=l 



OO. 



(43) 



Prom (l4*2j) . it is clear that X q > for all and g. In fact, if a X q = 0, then by the first complemen- 



tarity condition, p v q {k) = for all k, which is not possible by the last inequality in ([42 
Thus, it follows from the second complementarity condition that 

/ \ 



N 

E 

k=l 



log 



1 + 



\H qq (k)\ 2 p-(k) 



\ 



al(k) + Y,\H qr (k)\ 2 P^k) 



Rq + T v 



(44) 



which implies that the sequence {r u X q } is bounded for each q £ $7. We claim that lim J/ _ +D0 t u = 0. 
Otherwise, for some subsequence {t u : v £ k}, where k is an infinite index set, we have 
lim inf u r eK \_ too t u > 0. Thus, the subsequence {X q : v € k} is bounded for all q G The first 
complementarity condition in (|42l) then implies that {Pq(A:) : ^ E n} is bounded for all q € O and 
A; € A/". This is a contradiction to (l43j) . Therefore, the sequence {r^} J, 0. 

Consider now the normalized sequence {pV Hp"!!}, which must have at least one accumulation 
point; moreover, any such point must be nonzero. Let d°° be any such point. It is not difficult to show 
that d°° is a nonzero almost noiseless equilibrium. Moreover, from the inequality: 



JV 

E 

fc=i 



log 1 + 



\H qq (k)\ 2 p»(k) 



V 2 q(k)+E^q\Hqr(k)\* P » r (k) 



Rq < 0, 



(45) 



which is implied by dS|, it is equally easy to show that d°° £ 7 7 °(R ,tr ). Therefore, d°° is an element 
of A/"£o(R-*), which is a contradiction. This completes the proof of the existence of a solution to the 
problem (|32"T) . The boundedness of such solutions can be proved in a similar way via contradiction and 
by the same normalization argument. The details are not repeated. ■ 

Roughly speaking, the key condition Af£o(R*) = in the previous theorem is just the mathematical 
requirement that if the power p goes to infinity staying feasible, the system cannot approach a (noiseless) 
equilibrium. As such the previous theorem is rather natural, although it does not provide an effective 
way of checking the existence and boundedness of the solutions. To this end, however, we can now 
easily derive Theorem [5] from Theorem! 



A. 4 Proof of Theorem! 



In order to prove Theorem [5] we introduce a simple polyhedral set that will turn out to be a subset of 
P°"(R*). 

Poo(R*) 4 JJ { r (fc) e R« : Z fe (R*)r(fc) < o} , (46) 



fe=i 
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which is independent of the <x and where, we recall, the matrices Zfc(R*) are defined by (fTTI) . The key 
property for the existence Theorem [5] is stated in the following proposition. 

Proposition 24 Af£ (R*) C Poo(R*) n (M£ \ {0}). 

Proof. It suffices to note the following string of implications: 

^ |H M (fc)| 2 p,(A;) < (e^-l)^|tf gr (fc)| 2 Pr(fc), VgGO, (47) 

r+q 

o P e Poo(R*)- 

where the middle equivalence is by simple exponentiation. ■ 

It is known that, since each matrix Zfc(R*) is a Z-matrix, if each matrix Zfc(R*) is also a P-matrix, 
then a positive vector s(k) = (s q (k))^ =1 exists such that s T (/c)Zfc(R*) > [32j Theorem 3.3.4]. But 
this implies that Poo(R*) = {0}, and thus Af£o(R*) = 0. The first assertion of Theorem [5] then follows 
immediately from Theorem l23l It remains to establish the bound on the p(k) = (p q (k)) q=1 . This 
follows easily from the following two facts: 1) any solution p of (|32l) must belong to the set P°"(R*); 
2) a Z-matrix that is also a P-matrix must have a nonnegative inverse. ■ 



B Proof of Corollary [6] 

Consider the matrix Z max (R*) defined by (PT31) in Corollary [6j and assume it is a P-matrix. Set 
D(fe) = Diag {\H q q{k)\ 2 ^ =v The first assertion in the Corollary is immediate because 



Z fc (R*) > Z max (R*)D(/fc), VA: G M, 



(48) 



where the inequality is intended component-wise. Therefore, since all the matrices involved are Z- 
matrices, from the assumption on Z max (R*), it follows that all the Z^(R*) are also P-matricesQ Lets 
now prove the bounds (fl~5j) . Note first that 



/ af(k)(e R i -1) \ 



< 



maxa^(k) 



(49) 



1) J 



Furthermore we recall that, by (08]), we have [Z fc (R*)] 1 < [Z max (R*)D(fc)] 1 [32], and also that the 
inverse of a matrix that is P and Z is nonnegative [32j Theorem 3.11.10]. From all these facts, and 



7 The last statement can be easily proved using [32, Lemma 5.3.14]. 
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recalling (jT2j) , the following chain of inequalities easily follows for every k € N: 



D(fc) 



V Pg(k) J 



*\ l-i 



<D(fc)[Z fc (R*)] 



/ al(k){e R i-l) \ 



\ al(k)(e R h-l) ) 



< 



m&xa T .(k) 



which provides the desired bound (fT5"l) . 



[Z max (R*)] _1 (e R * - 1), 



(50) 



C Proof of Proposition [TT1 

We prove that the following three statements are equivalent for game when N = 1: 

(a) The game has a nonempty solution set; 

(b) The matrix Z(R*) is a P-matrix; 

(c) A^ (R*) = 0. 

(a) =4> (b). Since N = 1, any solution of (f32l) satisfies the equations 

/ \ 

\H qq (k)fp q (k) 



log 



1 + 



V 



al(k) + J2\HAk)\ 2 Pr(k) 

r^q 



which are easily seen to be equivalent to (|2l)|) . Since the right-hand constants in (|2T)|) are positive, it 
follows that a vector r* = (r*)^ =1 > 0, which is the solution of (|32"1) with N = 1, exists satisfying 
Z(R*)r* > 0. Hence (b) follows [32 Theorem 3.3.4], [HI Theorem 6.2.3]. The implications (b) => (c) 
^> (a) come directly from the more general case N > 1. Hence (a), (b), and (c) are equivalent. 

It remains to establish the assertion about the uniqueness of the solution. But this is clear because 
any one of the three statements (a), (b), or (c) implies that the matrix Z(R*) is a P-matrix, thus 
nonsingular [3T| Theorem 6.2.3]; and hence the system of linear equations (j26l) has a unique solution. 



D Proof of Theorem M 

The study of the uniqueness of the solution of the GNEP is complicated by the presence of a coupling 
among the feasible strategy sets of the users, due to the rate constraints. To overcome this difficulty 
we first introduce a change of variables of the game £f, from the power variables p q (k) to a set of rate 



22 



variables, in order to obtain an equivalent formulation of the original generalized Nash problem as a 
Variational Inequality (VI) problem, defined on the Cartesian product of the users' rate admissible 
sets. Then, building on this VI formulation, we derive sufficient conditions for the uniqueness of the 
GNE of the original game. It is important to remark that our VI formulation of the game <S differs 
from that of In fact, in [TT| the rate maximization game was formulated as an Affine VI defined 
on the Cartesian product of the users' power sets [TTJ Proposition 2]. Our VI instead, is defined by a 
nonlinear function, which significantly complicates the uniqueness analysis, as detailed next. 



D.l VI formulation 

Hereafter we assume that conditions of Theorem [5] are satisfied, so that a GNE of the game ?f is 
guaranteed to exist. 

Given the game 5f we introduce the following change of variables: 

\H qq (k)\ 2 p q (k) 



r q {k) 4 log 1 + 



<? 2 q (k)+Zr^\H q r(k)\ 2 Pr(k) 



, k G M, q G O, 



(51) 



with r q (k) satisfying the constraints 



N 



,(k) > 0, VA; G M, Vq G Q, and Yj r q (k) = R*, Vg G ft. 



(52) 



fc=i 



Observe that each r q (k) = if and only ifp q (k) = 0. Given r(k) = {r q (k))® =1 , let us define the Z-matrix 
Z k (r(k)) G R QxQ as: 



Z fc (r(fc)) 



\Hn(k)\ 2 -(e r ^ -l)\H 12 (k)\ 2 

-(e r ^-l)\H 21 (k)\ 2 \H 22 (k)\ 2 

.( e r Q (fc) _ i) \H Q1 (k)\ 2 -(e^« - 1) \H Q2 (k)\ 2 



-(e r ^-l)\H 1Q (k)\ 2 
-(e r ^-l)\H 2Q (k)\ 2 

\H QQ (k)\ 2 



(53) 

From (H2D, we have Z k (r(k)) > Z k (R*) for all k G J\f, where Z fc (R*) is defined in ((TTJ). It follows that 
each Z k (r(k)) is a P-matrix. 

According to (|5"T1) . the users' powers p(k) = (p q {k))® =1 are related to the rates r(k) by the following 
function 



/ p a (fc) \ 



\ PQ( k ) j 



( 



0(fc,r(fc))4 {Z^vikW 1 



a 2 (k)( e r ^ -1) 



{ a 2 Q (k){e^)-l) J 



fceJV. 



(54) 



Observe that (Z^(r(/c))) 1 is well-defined, since Z k (r(k)) is a P-matrix [32l Theorem 3.11.10]. 
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Using ([52]) and ([54"]) . and the fact that each p q (k) = if and only if r q (k) = 0, the KKT conditions 
of the Nash problem j6|) can be rewritten as (see (|32|) ): 



< r q (k) ± a 2 q {k) + J2 \H qr (k)\ 2 Mk,r(k)) - \H qq (k)\ 2 \ q > 0, \/k £ AT, Vq £ ft 

N r=1 (55) 



< X q , ^r q (k)=R* q , V ^ G0 



k=l 



where 4> r (k,r(k)) denotes the r-th component of cf>(k,r(k)), defined in (]5"4"]) . It is easy to see that ([55 
is equivalent to (note that as usual X q > for any solution of ([55])), VA; G TV, Vg G 



< r,(fc) 1 log a 2 (A;) + £ liT^AOPMMk)) - log (l# 3 #)| 2 ) + v q > 0, 

. V 7 (56) 



v q free, 



£r ff (fc) =^ 



fe=i 



Let us define 



* g (fc,r(*0) 4 log a 2 (A;) + £ \H qr (k)\ 2 r (fc, r(*)) - log (|if w (A;)| 2 ) , k £ N, q £ £1. 

V r=1 / (57) 

* ff (r)^(*,(fc,r(*)))£=i, *(r) 4 (<& g (r))J =1 , r4(r,)J =1 , 
with r 9 = (r q (k))f? =1 . Observe that each $> q (k,r(k)) in ([57]) is a well-defined continuously differen- 
tiable function on the rectangular box [0, R*] = Y\^=i [®iRq] c thus 3? in ([57]) is a well-defined 
continuously differentiable function on [0, R*]"^. 

Using ([571) . one can see that ([56]) are the KKT conditions of the VI (U, 3>) pU Proposition 1.3.4], 
where U is the Cartesian product of users' rate sets, defined as 



Q ( N ) 

U = ]J U q , where U q 4 ) Tq £ : £r g (*) = i?* L 
<?=i I fe=i J 



(5E 



and is the continuously differentiable function on [0,R*] defined in ([57]) . 

By definition, it follows that a tuple r* = (r*)® =l is a solution of the VI(f7, 3?) defined above if and 
only if, for all r q £ U q and q G fi, 



iV 



2 (r,(fc)-r*(fc); 



fc=i 



log a 2 (k) + g |^ r (A;)| 2 r (fc,rj£) - log {\H qq (k)\ 



r=l 



> 0. 



We rewrite now condition (1591) in a more useful form. To this end, let introduce 



T,(fc,p(fc)) 4 a 2 (k) + ^2\H qr (k)\ 2 Pr (k), k£N,q£tt, 

r+ q 



(59) 



(60) 



with p(fc) = (p 9 (A;))S =1 . For any solution p* = (p*(k))^_ 1 of ([3T]) (i.e., any GNE of the game ^ we 
have p* < p, where p = (p(k))^ =1 , with p(fc) = (p r (A0)^Li and i? r (/c) defined in (TT2]) . It follows that 

^ 2 (fc) < r,(A;,p*(fc)) < f g (fc) 4 <y 2 q (k) + Y^\H qr {k)\ 2 p r (k), VA; G jV, G O. (61) 
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Using (|60l) . we can write 



Q 



log\a 2 q (k) + Y,\H qr (k)\ 2 Pr(k) + \H qq (k)\ 2 Pq (k) ] = r q {k) + log (r q {k, p(fc)) ) , (62) 



so that condition (|59j) becomes 

A' 



E (r ff (fc)-r*(fc)) (log(r,(fc,p*(A;)))-log(|fr M (fc)| 2 )+r*(A;)) > 0, Vr, € Vg G 0, (63) 



fc=l 



where r q (k,p*(k)) is defined in ([60jl and p*(/c) = (p q {k)) q=1 , with each Pq(fc) = 4> q (k,r*(k)), and 
given in ([Ml) . Condition ([63]) will be instrumental for the study of the uniqueness of the GNE, as shown 
next. 



D.2 Uniqueness analysis 

Building on ([63]) . we derive now sufficient conditions for the uniqueness of the GNE of the game . 
Let pM 4 (gM)^, for i/ = 1,2, be any two solutions of the Nash problem in Given z/ = 1,2, 
k G A/", and g G H, let us define 



fJ")(A;) A lQg 1 + 



\H„{k)\*#\k) 



v 2 q (k) + Zr^\H qr (k)\ 2 &\k)J 

Adding the following two inequalities, which are obtained from the characterization ([63]) of a solution 
to the VI (£/,*): 

iV 

E (ff (k)-r^(k)) ( - log (|i/ w (fc)| 2 ) + log^Cfc)) + f «(A;) ) >0, Vqen, (65) 



and ^(^T ? (fc,pW(fc)). 



(64) 



and 



k=l 



N 



E fo (1) (fc)-ff (AO) (-log(|fl- w (A;)| 2 )+log(^ 2 ) (*0)+?f (A:)) > 0, G 0, (66) 



fe=i 

and rearranging terms, we obtain 

A^ 



N 



k=l 



k=l 



E W*) ^ E KfW-rW(A;) log(rW(A ; ))-log(T g ( 2 )(A ; )) 



< 



AT 



\ fe=i 



AT 



E ( l0 S(r q {1 \k)) - log(?f (£;))) Vq G O, 
\ fc=l 



which implies 



9(2) _ 9(1) 


< 


9 9 


2 



log(rf)-log(T«] 



Vg G n, 



(67) 



where r„^ = (Tq (k))^ =1 , for z/ = 1,2, and log(r^) has to be intended as the vector whose fc-th 
component is log(T q u \k)). 
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Invoking the mean- value theorem for the logarithmic function, we have that there exists some scalar 



s q (k) such that 
and, for each q € and k G A/", 



H qq (k)\ 



pf(fc) 



^ (2) (fc)-fW(A;) =log 1 + 



ff w (fc)| 2 p< 2) (fc) 



< s q {k) < \H qq (k)\< 



log 1 + 



r g {1) (k) 



(68) 



$\k) jff>(fc) \ |ff gg (fc)] 2 
T g (2) (£;) ? g (1) (fc) / 1 + 
$\k)-$\k)\ \H qq (k)\ 2 



ri 2 \k) 



+ 



\H qq (k)\ 2 $\k) 



^){k)-^\k)) 



-(1) /, s 
T q (k) rj L q 



\H qq {k)\ 2 



? q {2 \k)(l + s q (k)) 



Similarly, there exists some scalar uj q (k) such that 



f\k)<u q {k)<?^{k), 



(69) 



(70) 



and, for each q € and k G M, 

\og{T^{k))-\og(T q ^(k)) = i = -1— ^ jiy gr (fe)| 2 (p< 2 >(fc) -p< x >(fc) )- 



Introducing 



} {k) ±\H qq {k)\\&\k) -fc\k)), 



(71) 
(72) 



and using (|69j) and (|7T1) . the inequality ([67]) becomes 



E I I e,w + |g «SCf w E i^vwi 2 (sPw -rf"w 



fe=i 



r 9 ^ (^) r 



f 9 (2) (fc)(l+ S #)) 



(73) 



E (^E^)i 2 (^(fe)-^(*)) 

fe=l y « V ; r ^g 



By the triangle inequality and rearranging terms, from (|73p it follows 



N 

ilk U (2) (*0(i + *,(*oJ 



£ q (k) 



N / 

E TTTnEl^WI 2 (p^W-p^*) 



+ 



(74) 



N 

E 

fc=i 



|ff gg (fc)| 2 j4 1} (fc) 
r g (1) (fc) 



Y,\H q r(k)? {f 2 \k)-^\k) 
r+ q 



r q (2 \k)(l + s q (k)) 
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We bound now ([711) using the following: \/q £ Q, V/c € 

|^ ? (fc)| 2 #(fc) 



<e -Ji «-l, u = l,2, 



?<?\k) 
1 



l + s,(fc) 



1 1 1 

<^ <-^— , i/ = 1, 2, 



1 1 

< 



(75) 

(76) 
(77) 

(78) 



where (JTSJ) follows from $12D, §M from (EHI) and ([75]), ((77]) from JEE]), and J7E]) from ([701). Using 
GB-GB, (EH) can be bound as 



e -R* q 



iV 

E 



(79) 



< 



N / R* -, 

i e i — 1 



\fc=l \ 9 V ^ r^q 



+ 



e(^Ei^w(# 

fc=l y <? v ; 



(80) 



e(^eiw(#' 

fc=l y 1 K ' r^q 



P?\k)-p£\k) 



e R* q 



fc=l Vr^g 



£r(fc) 
f r (fc) 



— E ( max /V(fc) 



AT / /, n \ 2 



(82) 



33) 



where: §79$) follows from (176 )1 and (1771): (|8f) |) follows from (JZU), ((771) and {ZSJ); and in (|82| we have 
defined 

M*0 = ^j^^jL 9,r e ft, and k <EJ\T. (84) 



Let B = B(R*) = [b 9 r]^ r=1 be the nonnegative matrix, where 



'qr 



if i = j 



e i max P„ r (k) if i / j, 
fceA/" 



(85) 



and let B be the "comparison matrix" of B, i.e., the matrix whose diagonal entries are the same as 
those of B and the off-diagonal entries are the negatives of those of B (see (fl7l) ). Note that B is a 
Z-matrix. 
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Introducing 




with £q(k) denned in (|72j) . and concatenating the inequalities in (|79j) for all q E fi, we deduce 

Bt < 0, (87) 

with t = (t q )® =1 . If B is a P-matrix, then it must have a nonnegative inverse [32J Theorem 3.11.10]. 
Thus, by (15Tj) . we have t < 0, which yields t = 0. This proves the uniqueness of the GNE, under 
conditions of Theorem [9] 

Remark 25 An alternative approach to establish the uniqueness of the solution of the Nash problem 
([!]) is to show that under a similar hypothesis, the function <&(r) in (f571) is a "uniformly P-function" on 
the Cartesian product set U. In turn, the latter can be proved by showing that the Jacobian matrix 
J3>(r) of the function <&(r) is a "partitioned P-matrix" uniformly for all r E U. We adopt the above 
proof because it can be used directly in the convergence analysis of the distributed algorithm to be 
presented subsequently. 



E Proof of Corollary EE 



To prove the desired sufficient conditions for B(R*) in (fT7]) to be a P-matrix, we use the bounds ([15 
in Corollary [6] as shown next. 

We provide first an upper bound of each f3 qr (k), defined in ([M]) . Let d = (d q )® =1 be defined as 

d 4 (Z max (R*))" 1 (e R * - l) , (88 



with Z max (R*) given in $T5§ and R* = (R* q )® =1 . By (EJ, we have 



f q (k) < a 2 Jk)+( maxa 2 ,(fc)) y \Esi^d T , 
" q Vr'en rK V f^ q \H rr (k)\ 2 



\/k E Af, Vq, red, and q / 



(89) 



where r q {k) is defined in (|6T]l . Introducing (|89l) in (|84|) . we obtain 



|ggr(fe)| 

| | 



of (fc) + V™ of(*) ) fy q 



\H qr '{k)\ - 

TU, r l 



\H r / r /(k)\ 



2 



VA: E jV, Vg, r E f2, and q ^ r. 

(90) 

Hence, recalling the definitions of the constants x (assumed in (0, 1)) and p > 1, as given in (fT9j) and 
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([20]) . respectively, we deduce 

n „n / l^r(fc)| 2, 

max. Bar (k) < max — - ^ 

■ ■ " ' ~ 1 keN \H rr (k)\\ 



max 



+ 



( 

max 



max 



r V? 



|-£V'(fc)l' 

max J 

kcAf \H r > r >(k)\ 



<- flmax 
— 



a 2 r (k)\ 
max | + 



max max — „ . 
r'esj VfceA" cr|(fc) 



X 



1 J > , Vg, r£!l, and q ^ r, 



(91) 

where /3™ ax is denned in (fT8l) . 

Prom ([9T1) . one infers that condition (j23j) implies that B(R*) in (jXTJ) is diagonally dominant, which 
is sufficient to guarantee the P-property of B(R*) [HH Theorem 6.2.3], and thus the uniqueness of the 
GNP of © (Theorem ED. 

It remains to show that (|2~1) is equivalent to (|23l) if (Tg(k) = a^(k), \/r,q € Q and k G jV. In this 
case, (|23l) reduces to 

PE^ X < Xe" 2fl ', VgGO, (92) 



or equivalently 



/miax 
gen 



(e^-l)^/?, 



max 
qr 



rj^q 



< x max ( e i " 9 — { """ 
q en 



(93) 



which is clearly equivalent to J2 



F Proof of Theorem [14] and Theorem [18 

The proof of the convergence of both the sequential and simultaneous IWFAs is similar to the proof of 
the uniqueness of the GNE of the game £f as given in Appendix [Dj The difference is that instead of 
working with two solutions of the GNEP (and showing that they are equal under certain conditions), 
we consider the users' power allocation vectors produced by the algorithms in two consecutive iterations 
and derive conditions under which their respective distances to the unique solution of the game contract 
under some norm. 

We focus first on the convergence of the simultaneous IWFA. Then, we briefly show that a similar 
analysis can be carried out also for the sequential IWFA. Throughout the following analysis, we assume 
that conditions of Theorem [9] are satisfied. 
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F.l Convergence of simultaneous IWFA 

Let us define p( n+1 ) = (p q n ) q=1 , with p q n+1 ^ = {Pq n+1 \k))^ =l denoting the power allocation vector 
of user q at iteration n + 1 of the simultaneous IWFA given in Algorithm 1, and let 

|2 Jn+i), 



„(n+l) 



(k) ± log 1 + 



1 \H qq {k)\ 2 P y r L '(k) 



G Af, gell, n = 0,1,... , 



(94) 



with 4 n \k) = T q (k,p^(k)), where p (n 1(k) = {p q n \k))® =1 and r 5 (fc, p (n) (*0) is defined as in j60 



According to the simultaneous IWFA, at iteration n + 1, the power allocation Pg™ +1 ^ of each user q 
must satisfy the single-user waterfilling solution J7)) (see (|30T) ). given the allocations pi, = (Pr n ^)^ g= iof 
the other users at the previous iteration. It follows that each p q n+1 ^ and r q n+l \k) satisfy (see ([56]) ) 



< r q n+l \k) ± log(r 9 w (fc)) + r^'ik) - log(\H qq (k)\ 2 ) + u q n ^> > 0, k G Af, 



>+l) 



(n+1) 



4 n+1) free, ^ r> +1) (£;) = i?; 



fc=l 



Vn = 0, 1, . . . , 
(95) 



or equivalently (see ([M]) ). for all c/ G Q and n = 0, 1, 



E [r q (k)-ri n+1) (k)) (-log(\H qq (k)\ 2 )+log(T^(k))+r q n+1 \k)) > 0, Vr g G 17, , (96) 
fe=i 

where f7 g is defined in (|55j) . Let p* = (p^)^ =1 denote the unique GNE of the game <S (i.e., the unique 
solution of ([32]) ) and r* = (r*)^ =1 the unique rate solution of (|55l) . Note that each r* satisfies the 
constraints in j96j). Hence, following the same approach as in Appendix [D] to obtain (|67|) from ([63]) . we 
deduce 



r ("+l) _ r * 


< 


1 9 


2 



log(rW)-log(r*) , Vgefl, Vn = 0,l,..., 



(97) 



where t* = (T q (k))^ =l and T q {k) = r g (/c, p*(k)). Similarly to ([69]) . there exists some scalar s q n \k) such 
that 



and, for all G A/", g G fi, and n = 0, 1, . . . 



T-Jk) 



((n+1) / , n 
(pt +1 \k) -p* q (k)) +^y-^Ei^( fc )l 2 - #(*) 



|77 gg (fc)| 2 



(99) 



r*(k)(l + s q n \k)) 



Moreover, there exists some scalar uj q n \k) such that 



f\k)<u^{k)<rUk), 



(100) 
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and, for all k € A/", q £ ft, n = 0, 1, . . . , 

Iog(rW(fc))-log(r;(fc)) = -4— U n Hk)-r*(k)) = —^J2\ H ^)\ 2 U n) (k)-p;(kf 

J (k) uj^ >{k) T+q 

(101) 

Introducing the vector t( n ) = (t q n ^)® =l , with t q defined as 



An) A 



2 



^ '-J. 1^,n = ,l , (102) 

with f q (k) given in fHH) and using (ESJ) and (fTUil . (j9TD leads to 

(DiaglBCR^)}) t( n+1 ) < (off-Diag{B(R' t )})t( n ), n = 0, 1, . . . ., (103) 

where Diag{B(R*)} and off-Diag {B(R*)} are the diagonal and the off-diagonal parts of B(R*), re- 
spectively, with B(R*) defined in (j85l) . Based on (|103p . the proof of convergence of the simultaneous 
IWFA is guaranteed under conditions on Theorem (H as argued next. 

According to (32J, Lemma 5.3.14], the P-property of B(R*), with B(R*) defined in ([T7j) . is equivalent 
to the spectral condition: 

p [(Diag{B(R*)}) - 1 off-Diag{B(R*)}] < 1, (104) 

where p(A) denotes the spectral radius of A. Therefore, by (|103|) and |104p . the sequence {t( n )} 
contracts under a certain matrix norm; hence it converges to zero. The claimed convergence of the 
sequence {p^} follows readily. ■ 

F.2 Convergence of sequential IWFA 



The convergence of the sequential IWFA described in Algorithm 2 can be studied using the same 
approach as for the simultaneous IWFA. The difference is the final relationship between the error vectors 
in two consecutive iterations of the algorithm. More specifically, using the vectors = (t q n) )® =1 , with 
t q defined as in (|102p . one can see that the sequential IWFA leads to 



(Diag{B(R*)}-Low{B(R*)})1> +1 ) < (Up{B(R*)}) t^ n \ n = 0,l,...., (105) 

where Low{B(R*)} and Up{B(R*)} denotes the strictly lower and strictly upper triangular parts of 
B(R*), respectively. The above inequality implies 

t ("+i) < (Diag{B(R*)}-Low{B(R*)})- 1 (Up {B(R*)} ) t< n > = Tt^, n = 0,l,...., (106) 

where 

T 4 (Diag {B(R*)} - Low {B(R*)} ) _1 Up {B(R*)} . (107) 
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In (|105l) we used the fact that, under the P-property of the Z-matrix Diag {B(R*)} — Low {B(R*)} (due 
to the fact that all its principal minors are less than one), the inverse (Diag {B(R*)} — Low {B(R*)} ) _1 
is well-defined and nonnegative entry-wise [321 Theorem 3.11.10]. 

According to (|106p . the convergence of the sequential IWFA is guaranteed under the spectral con- 
dition 

p(T)<l, (108) 
which is equivalent to the P-property of B(R*) [32J Lemma 5.3.14], with B(R*) defined in (fl7l) . ■ 
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BS 



(a) 




(b) 

Figure 1: Probability of existence (red line curves) and uniqueness (blue line curves) of the GNE versus d 
[subplot (b)] for a 7-cell (downlink) cellular system [subplot (a)] and rate profiles R* q — 1 bit/symb/subchannel 
(square markers) and R* = 2 bit/symb/subchannel (cross markers), Vq S fi. 



35 




iterations (n) 



(a) 




iterations (n) 



(b) 

Figure 2: Rates of the users versus iterations: sequential IWFA (solid line curves), simultaneous IWFA (dashed 



line curves), Q = 10, d rq /d qr , d rr 



1, 7= 2.5. 
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